arXiv: 1509.0208Ivl [astro-ph.CO] 7 Sep 2015 


The Rethermalizing Bose-Einstein Condensate of 
Dark Matter Axions 

Nilanjan Banik, Adam Christopherson, Pierre Sikivie, Elisa Maria Todarello 
University of Florida, Gainesville, FL 32611, USA 

DOI: will be assigned 

The axions produced during the QCD phase transition by vacuum realignment, string 
decay and domain wall decay thermalize as a result of their gravitational self-interactions 
when the photon temperature is approximately 500 eV. They then form a Bose-Einstein 
condensate (BEC). Because the axion BEG rethermalizes on time scales shorter than the 
age of the universe, it has properties that distinguish it from other forms of cold dark 
matter. The observational evidence for caustic rings of dark matter in galactic halos is 
explained if the dark matter is axions, at least in part, but not if the dark matter is entirely 
WlMPs or sterile neutrinos. 


1 Axion dark matter 


The story we tell applies to any scalar or pseudo-scalar dark matter produced in the early 
universe by vacuum realignment and/or the related processes of string and domain wall decay. 
However, the best motivated particle with those properties is the QCD axion since it is not only 
a cold dark matter candidate but also solves the strong CP problem of the standard model of 
elementary particles mm- So, for the sake of definiteness, we discuss the specific case of the 
QCD axion. 

The Lagrangian density for the axion field 4>{x) may be written as 

+ - ( 1 ) 


where the dots represent interactions of the axion with the known particles. The properties 
of the axion are mainly determined by one parameter / with dimension of energy, called the 
‘axion decay constant’. In particular the axion mass is 
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in terms of the pion decay constant /,r, the pion mass and the masses m„ and of the 
up and down quarks, and the axion self-coupling is 
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All couplings of the axion are inversely proportional to /. When the axion was first proposed, / 
was thought to be of order the electroweak scale, but its value is in fact arbitrary [5]. However 
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the combined limits from unsuccessful searches for the axion in particle and nuclear physics 
experiments and from stellar evolution imply / > 3 • 10® GeV [4]. 

An upper limit / < 10^® GeV is obtained from the requirement that axions are not 
overproduced in the early universe by the vacuum realignment mechanism [^, which may 
be briefly described as follows. The non-perturbative QGD effects that give the axion its 
mass turn on at a temperature of order 1 GeV. The critical time, defined by m(ti)ti = 1, is 
~ 2 • 10“^ sec(//10^® GeV)3. Before ti, the axion field ^ has magnitude of order /. After 
ti, (j) oscillates with decreasing amplitude, consistent with axion number conservation. The 
number density of axions produced by vacuum realignment is 
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where a{t) is the cosmological scale factor. Their contribution to the energy density today 
equals the observed density of cold dark matter when the axion mass is of order 10“® eV, with 
large uncertainties. The axions produced by vacuum realignment are a form of cold dark matter 
because they are non-relativistic soon after their production at time G. Indeed their typical 
momenta at time ti are of order l/G, and vary as l/a(t), so that their velocity dispersion is 
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The average quantum state occupation number of the cold axions is therefore 
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Af is time-independent, in agreement with Liouville’s theorem. Gonsidering that the axions are 
highly degenerate, it is natural to ask whether they form a Bose-Einstein condensate nm. We 
discuss the process of Bose-Einstein condensation and its implications in the next section. 

The thermalization and Bose-Einstein condensation of cold dark matter axions is also dis¬ 
cussed in refs. [iliiiniin] with conclusions that do not necessarily coincide with ours in all 
respects. 


2 Bose-Einstein condensation 

Bose-Einstein condensation occurs in a fluid made up of a huge number of particles if four con¬ 
ditions are satisfied: 1) the particles are identical bosons, 2) their number is conserved, 3) they 
are highly degenerate, i.e. Af is much larger than one, and 4) they are in thermal equilibrium. 
Axion number is effectively conserved because all axion number changing processes, such as 
axion decay to two photons, occur on time scales vastly longer than the age of the universe. 
So the axions produced by vacuum realignment clearly satisfy the first three conditions. The 
fourth condition is not obviously satisfied since the axion is very weakly coupled. In contrast, 
for Bose-Einstein condensation in atoms, the fourth condition is readily satisfied whereas the 
third is hard to achieve. The fourth condition is a matter of time scales. Gonsider a fluid that 
satisfies the first three conditions and has a finite, albeit perhaps very long, thermal relaxation 
time scale r. Then, on time scales short compared to r and length scales large compared to a 
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certain Jeans’ length (see below) the fluid behaves like cold dark matter (CDM), but on time 
scales large compared to r, the fluid behaves differently from CDM. 

Indeed, on time scales short compared to r, the fluid behaves as a classical scalar field since 
it is highly degenerate. In the non-relativistic limit, appropriate for axions, a classical scalar 
field is mapped onto a wavefunction ip by 

(^(r, t) = ■ (7) 

The field equation for ^(x) implies the Schrodinger-Gross- Pitaevskii equation for ip 

idtip =+ V{r,t)ip (8) 

2m 

where the potential energy is determined by the fluid itself: 

V{r,t) =m<^{r,t)--^\ilj{r,t)\‘^ . (9) 

The first term is due to the fluid’s gravitational self-interactions. The gravitational potential 
<I>(r, f) solves the Poisson equation: 


= AirGinn 


( 10 ) 


where n = |'0p. The fluid described by ip has density n and velocity ii’ = &vg{ip). Eq. ([5]) 
implies that n and v satisfy the continuity equation and the Euler-like equation 

dtv + {v ■ W)v =——VV — Vq (11) 
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where 

1 

^ 2m^ y/n 
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q is commonly referred to as ‘quantum pressure’. The Vg term in Eq. dlB is a consequence of 
the Heisenberg uncertainty principle and accounts, for example, for the intrinsic tendency of a 
wavepacket to spread. It implies a Jeans length [T2] 
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where p = nm is the energy density. On distance scales large compared to £j, quantum pressure 
is negligible. CDM satisfies the continuity equation, the Poisson equation, and Eq. (ITT]) without 
the quantum pressure term. So, on distance scales large compared to £j and time scales short 
compared to r, a degenerate non-relativistic fluid of bosons satisfies the same equations as CDM 
and hence behaves as CDM. The wavefunction describing density perturbations in the linear 
regime is given in ref. m- 

On time scales large compared to r, the fluid of degenerate bosons does not behave like 
CDM since it thermalizes and forms a BEC. Most of the particles go to the lowest energy state 
available to them through their thermalizing interactions. This behavior is not described by 
classical field theory and is different from that of CDM. When thermalizing, classical fields 
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suffer from a ultraviolet catastrophe because the state of highest entropy is one in which each 
held mode has average energy ksT where T is temperature. In contrast, for the quantum held, 
the average energy of each mode is given by the Bose-Einstein distribution, and the ultraviolet 
catastrophe is removed. To see whether Bose-Einstein condensation is relevant to axions one 
must estimate the relaxation rate E = i of the axion huid. We do this in the next section. 

When the mass is of order 10“^^ eV and smaller, the Jeans length is long enough to affect 
structure formation in an observable way M- Because we are focussed on the properties of 
QCD axions, we do not consider this interesting possibility here. 


3 Thermalization rate 


It is convenient to introduce a cubic box of size L with periodic boundary conditions. In the 
non-relativistic limit, the Hamiltonian for the axion huid in such a box has the form 

H = ^^^ujjajaj + 'y ^ (14) 

j 

with the oscillator label j being the allowed particle momenta in the box p = ^ {rix,ny,nz), 
with rix, Uy and integers, and the given by [7] 

\P3,P4, _ \ P3,Pi A P3,Pi 

pi,p2 Pl,P2 

where the hrst term 

A P3,Pi _ ^ A- - - - 

P1.P2 “ P1+P2.P3-I-P4 

is due to the self-interactions, and the second term 


P4P 

is due to the gravitational self-interactions. 

In the particle kinetic regime, dehned by the condition that the relaxation rate F = ^ is 
small compared to the energy dispersion 6uj of the oscillators, the Hamiltonian of Eq. (1141) 
implies the evolution equation 
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Ml — 2 \MiMj{Mi + 1)(A4 + 1) ~ MiMkiMi + 1)(A/) -l- 1)] 2'K5{uJi + uij — u>k — W;) 

(18) 

for the quantum state occupation number operators A/) (t) = aj The thermalization rate 

in the particle kinetic regime, is obtained by carrying out the sums in Eq. (IlSp and estimating 
the time scale r over which the Mj change completely. This yields [6l [7] 


T ^ n (T 5v M 


(19) 


where cr is the scattering cross-section associated with the interaction, and M is the average state 
occupation number of those states that are highly occupied. The cross-section for scattering 
by \(jr self-interactions is a\ = ■ For gravitational self-interactions, one must take the 
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cross-section for large angle scattering only, ag ~ i since forward scattering does not 

change the momentum distribution. 

However, the axion fluid does not thermalize in the particle kinetic regime. It thermalizes in 
the opposite “condensed regime” defined by F >> 6uj. In the condensed regime, the relaxation 
rate due to Xcj)'^ self-interactions is [SE] 


r 


nX 

Am? 


and that due to gravitational self-interactions is 


Fg ~ iirGnm^i^ 


( 20 ) 


( 21 ) 


where i = is, as before, the correlation length of the particles. One can show that the 
expressions for the relaxation rates in the condensed regime agree with those in the particle 
kinetic regime at the boundary i5a; = F. 

We apply Eqs. dini) and dSD to the fluid of cold dark matter axions described at the end 
of Section 1. One finds that Fa)^) becomes of order the Hubble rate, and therefore the axions 
briefly thermalize as a result of their X(j)^ interactions, immediately after they are produced 
during the QCD phase transition. This brief period of thermalization has no known impli¬ 
cations for observation. However, the axion fluid thermalizes again due to its gravitational 
self-interactions when the photon temperature is approximately [iE] 


Fbec 


500 eV 


/ 
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( 22 ) 


The axion fluid forms a BEG then. After BEG formation, the correlation length £ increases till 
it is of order the horizon and thermalization occurs on ever shorter time scales relative to the 
age of the universe. 


4 Observational consequences 

As was emphasized in Section 3, the axion fluid behaves differently from CDM when it ther¬ 
malizes. Indeed when all four conditions for Bose-Einstein condensation the axions almost all 
go to their lowest energy available state. CDM does not do that. One can readily show that in 
first order of perturbation theory and within the horizon the axion fluid does not rethermalize 
and hence behaves like CDM. This is important because the cosmic microwave background 
observations provide very strong constraints in this arena and the constraints are consistent 
with CDM. In second order of perturbation theory and higher, axions generally behave differ¬ 
ently from CDM. The rethermalization of the axion BEC is sufficiently fast that axions that 
are about to fall into a galactic gravitational potential well go to their lowest energy state 
consistent with the total angular momentum they acquired from nearby protogalaxies through 
tidal torquing [7]. That state is a state of net overall rotation. In contrast, CDM falls into 
galactic gravitational potential wells with an irrotational velocity field. The inner caustics are 
different in the two cases. In the case of net overall rotation, the inner caustics are rings m 
whose cross-section is a section of the elliptic umbilic D _4 catastrophe m, called caustic rings 
for short. If the velocity field of the infalling particles is irrotational, the inner caustics have 
a ‘tent-like’ structure which is described in detail in ref. m and which is quite distinct from 
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caustic rings. There is observational evidence for caustic rings [18]. It was shown [19] that the 
assumption that the dark matter is axions explains not only the existence of caustic rings but 
also their detailed properties, in particular the pattern of caustic ring radii and their overall size. 
Furthermore, it was shown [20] that axion BEC solves the galactic angular momentum problem, 
the tendency of CDM to produce halos that are too concentrated at the center compared to 
observations. 

In a recent paper m. j- Dumas et al. compare the predictions of the caustic ring model 
with the rotation curve of the Milky Way and the observations of the Sagittarius sattelite’s 
tidal disruption. 
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